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Abstract — In this paper we describe the effect of imperfections 
in the radio environment map (REM) information on the per- 
formance of cognitive radio (CR) systems. Via simulations we 
explore the relationship between the required precision of the 
REM and various channel/system properties. For example, the 
degree of spatial correlation in the shadow fading is a key factor 
as is the interference constraint employed by the primary user. 
Based on the CR interferers obtained from the simulations, we 
characterize the temporal behavior of such systems by computing 
the level crossing rates (LCRs) of the cumulative interference 
represented by these CRs. This evaluates the effect of short term 
fluctuations above acceptable interference levels due to the fast 
fading. We derive analytical formulae for the LCRs in Rayleigh 
and Rician fast fading conditions. The analytical results are 
verified by Monte Carlo simulations. 

I. Introduction 

It is now well known [1], [2] that granting exclusive licences 
to service providers for particular frequency bands has led to 
severe under-utilization of the radio frequency (RF) spectrum. 
This has led to global interest in the concept of cognitive radios 
(CRs) or secondary users (SUs). These CRs are deemed to 
be intelligent agents capable of making opportunistic use of 
radio spectrum while simultaneously existing with the legacy 
primary users (PUs) without harming their operation. 

In addition to ensuring quality of service (QoS) operation, 
the most important and challenging task for the CRs is to 
avoid adverse interference to the incumbent PUs. Hence, it is 
necessary to develop schemes that can help PUs avoid such 
harmful interference. In addition to the primary exclusion zone 
(PEZ) [3] approach, the recently developed [4] methods based 
on radio environment maps (REMs) [5], [6] can also help 
achieve this goal. In an earlier paper [4] we have shown that 
under certain conditions the REM based approach can result in 
substantially higher numbers of permissible CRs than the PEZ 
approach. Hence, the focus of this paper is an REM scheme. 
In particular we consider a REM which stores signal strength 
data from point to point in a regular grid. CRs have access to 
this REM and can therefore evaluate their impact on the PU 
and maintain acceptable interference levels as long as they can 
obtain positional information on the other CRs and the PU. 

The REM based approaches heavily depend on "quantity" 
and "quality" of the REM information available. Defects in 
REM information can seriously affect the PU performance. In 
a similar manner, temporal variations in the CRs' interfering 
signals can degrade the PU performance even though the CR 



level may be acceptable on average. These two aspects form 
the focus of this paper. In both situations we assume that the 
PU is willing to suffer some reduction in SNR, so that an 
allowable level of interference is provided to enable the CR 
operation. In particular, we make the following contributions: 

• We determine the impact of coarse REM information. 
We show that when the REM for a given area is dis- 
cretized then the total CR interference is significantly 
underestimated when realistic grid sizes are considered. 
For example, for a grid size of 50 m x 50 m, the actual 
SINR is worse than the target SINR by at least 1 dB 
for 8% of the time. We also determine the interaction 
between shadow fading correlation and REM grid size 
and evaluate their impact on interference estimation. 

• We determine the level crossing rate (LCR) and ex- 
ceedance duration (AED) of the CR-PU interference for 
a number of scenarios including Rayleigh fading, Rician 
fading, and various CR interferer profiles. The LCR is 
determined via analysis and confirmed via simulation. 
Results show that the LCR is maximum at or around the 
maximum interference threshold and is virtually zero 5 
dB beyond this point. We also show that for urban areas 
which are characterized by a strong LOS component, the 
interference rarely crosses the threshold and when it does, 
it only exceeds the threshold value for small duration. 

The rest of the paper is organized as follows: Section II 
describes the system model and the REM. Section III char- 
acterizes the instantaneous composite CR interference to the 
PU system in terms of the LCRs. In Section IV we present 
simulation and analytical results. Finally, in Section V we 
describe our conclusions. 

II. System Model and REM 

Consider a PU receiver in the center of a circular region 
of radius R. The PU transmitter is located uniformly in an 
annulus of outer radius R and inner radius Rq centered on the 
PU receiver. It is to be noted that we place the PU receiver 
at the center only for the sake of mathematical convenience. 
The use of the annulus restricts devices from being too close 
to the receiver. This matches physical reality and also avoids 
problems with the classical inverse power law relationship 
between signal strength and distance [7]. In particular, having 
a minimum distance, R , prevents the signal strength from 
becoming infinite as the transmitter approaches the receiver. 



Similarly, we assume that multiple CR transmitters are uni- 
formly located in the annulus. At any given time, each CR 
has a probability of seeking a connection, given by the activity 
factor, p. The number of CRs wishing to operate is denoted 
Ncr- Of these CRs, a certain number will be accepted 
depending on the allocation mechanism. Hence, a random 
number of CRs denoted N < Ncn will transmit during the 
PU transmission and create interference at the PU receiver. 

The received signal strength for both the PU transmitter to 
PU receiver and CR transmitter to PU receiver is assumed to 
follow the classical distance dependent, lognormal shadowing 
model. For a generic interferer, this is given by 



I = BLr~~< = BW x/w r-~< = 



Be x r~ 



(1) 



where r is the random distance from the transmitter to the 
receiver, 7 is the path loss exponent (normally in the range 
of 2 to 4) and L is a shadow fading variable. The lognormal 
variable, L, is given in terms of the zero mean Gaussian, X, 
which has standard deviation a (dB) and X = j3X where 
(3 = log(10)/10. The standard deviation of X is denoted by 
a x . The constant B is determined by the transmit power. The 
desired primary signal strength, S, has the same form, with a 
different transmit power, so that S = AL p r~^ . The constant 
A is determined so as to give an SNR greater than 5 dB, 95% 
of the time. However, the constant B depends upon A, 7 and 
the ratio of R and R c as given in [4] . Note that all the links are 
assumed to be independent and identically distributed (i.i.d.) 
so that spatial correlation is ignored. 

A. A Perfect REM 

A REM can hold a wide variety of information [6] and 
it is not clearly understood at present what constitutes a 
practical and effective REM. In this work we assume that 
the REM contains signal strength data. In a perfect REM the 
signal strength from all source coordinates to all destination 
coordinates is known. With this perfect REM a CR controller 
[4] can select those CRs for operation which satisfy a given 
interference constraint. The CR controller requires positional 
information for the PU and the CRs, and can then use the 
REM to compute the overall SINR of the PU where 

SINR = (2) 

In (O, S is the signal strength of the PU, a 2 is the noise 
power and i s tne a gg re 8 ate interference of the N 

selected CRs. The interference constraint used is that the CR 
interference must not reduce the PU SNR by more than 2 dB. 
All results shown in the paper are for 2 dB buffer. The value 
of 2 dB was chosen arbitrarily and an exploration of effects 
of this buffer will appear in a future work. In this paper the 
centralized approach in [4] is used to select the N CRs for 
operation. 

B. Modeling of REM Imperfections 

In practise a perfect REM is impossible and for practical 
purposes the REM information is discretized in the form of 
grid points with grid size, A. Hence, the central controller 



allocating CRs will formulate its decisions on the basis of 
REM information obtained from the grid points, rather than 
from exact signal strength data. Hence, an interfering signal 
strength, /, will be estimated by / from the REM. The estimate 
is obtained from the grid-to-grid path in the REM which is 
closest to the actual signal path. 

We consider the CR signal strength to be of the form given 
in ([T). The REM predicted signal strength is given by: 



Be x f-~< 



(3) 



where f is the distance between the transmitter and the receiver 
in the REM grid and X is correlated with X by: 



X = pX + s/l - p 2 E 



(4) 



In (O E is independent and identically distributed (i.i.d.) 
with X. Assuming a distance, di, between the actual and 
REM based position of the CR and a distance, d p , between 
the actual and REM based location of the PU receiver, the 
correlation coefficient p can be obtained using an extension of 
Gudmundson's model [8] as: 



p = 0.h d * /Dd x 0.5 dp/Dd 



(5) 



In ([5]) Dd is the so called decorrelation distance i.e., the 
distance at which the correlation between X and X drops 
to 0.5. The effect of flawed REM information on the signal 
strength between the primary transmitter and its receiver can 
also be modeled using ([3]), (01 and (0. Simulation results 
of this model based on parameter values of a suburban 
macrocellular environment are given in Section IV. 

III. Instantaneous CR performance 

The CR allocation policy is based on mean signal and 
interference levels. As a result, even if the 2 dB buffer 
is exactly met the instantaneous fast fading will result in 
fluctuations of the SINR both above and below the buffer. It is 
therefore of interest to investigate how often and how long the 
SINR exceeds the buffer. As a first look at this problem we fix 
the PU signal power and consider the instantaneous variation 
of the interference only. In this scenario the 2 dB SINR 
buffer becomes a threshold of -2.33 dB for the interference 
(as shown in Figs. 4-7). Hence in this section we focus on the 
instantaneous temporal behavior of the aggregate interference. 
For this purpose we evaluate the LCR (and thus the average 
exceedance duration (AED)) of the cumulative interference 
offered by the CRs obtained using the centralized approach 
of [4] with imperfect REM information. First we calculate 
the LCRs for Rayleigh environment and then we characterize 
them for Rician fading conditions. In future work the full 
temporal behavior of the SINR should be considered, but this 
preliminary investigation still yields useful results and insights. 

A. LCRs for Rayleigh Fading 

For a given set of CR interferers, the instantaneous aggre- 
gate interference under Rayleigh fading, lR ay {i), is given by: 



N 



I Ray {t)=J2 I *M t )\' 



(6) 



where li represents the interference power of the i th CR, 
|/ii(i)| 2 is a standard exponential random variable with unit 
mean and N is the number of interfering CRs. From (O, 
the aggregate interference is represented as a weighted sum 
of exponential random variables. Such weighted sums can 
be approximated by a gamma variable [9]. Simulated results 
show that the gamma fit is very good, but are not shown here 
for reasons of space. It should be noted that the exact LCR 
computation for such weighted sums was given in [10] for 
the case of three and four branch maximal ratio combining 
(MRC) by providing special function integrals. Recently, more 
general expressions for arbitrary number of branches have 
been derived in [11]. However, the approach of [ 1 1 ] results in 
numerical difficulties, especially for large values of N, which 
can be the case for CR systems. Hence an approximation is 
useful to overcome these problems and to provide a much 
simpler solution. A gamma variable with shape parameter r 
and scale parameter 9 has a mean and variance given by r/9 
and r/9 2 respectively and probability density function (PDF): 



f(x) = T{r)- 1 9 r x r - 1 exp(-9x), 



x>0 



(7) 



Thus, approximate LCRs for © can be found by calculating 
the LCR of the equivalent gamma process. The LCR for a 
gamma process has been calculated in [12]. Thus, the crossing 
rate of I Ray (t) across a threshold, T, can be approximated by: 



LCR /HoB (T) 



2T(r) V 7T 

E(I Ra y{t)) 2 /Var(I Ray (t)), 



where r = 

E(I Ray (t))/Var(I Ray (t)) and R(0) = p Ray (0) is the 
second derivative of the autocorrelation function (ACF) of 
I Ray (t). Hence, to compute the LCR in (0 only the mean, 
variance and ACF of the random process in (0 are required. 

The first two moments of (O are simple to compute as 



E(lRay(t)) = Z? =1 Ii and Var(I Ray (t)) = Eti 1 !- To 
calculate the ACF, note that hi(t + t) can be written as: 



hi(t + T )= Pl {T)h l {t) + y/(l-fi(T))ei(t), (9) 

where ei(t) is independent of hi{i) and statistically identical 
to hi(t). Assuming a Jakes' fading process, pi(r) is the zeroth 
order Bessel function of the first kind, Jo(27t/dt) and jo is 
the Doppler frequency. Using © we have: 

N 

E[I Ray (t)I Ray (t + r)} = ]T LJ 1 E[\Ht)\ 2 \h^ + ^)\ 2 ] 

N , N 

i=£j V i=\ 

xM*)l 2 + (i-^(r))Mt)| 2 )] 

TV N N 

= E^+E J *+E^) 

ijkj i—1 i—1 

/ N ^2 N 

= (E J +E J ^)> ("» 



i=l 



where in the second to last step above, we have used the fact 
that cross products have zero mean and that E^/i^i)! 4 ] = 2. 
The ACF of (|6]l is given by: 

E(I Ray {t)I Ray (t + T))~E(I Ray (t))E(I Ray (t + r)) 



^/Var{I Ray {t))Var{I Ray {t + r)) 



and with the relevant substitutions, the ACF becomes: 



PRay{T) 



T N I 2 



(11) 



(12) 



Finally, using the expansion J (2-Kf^T) = 1 — 7r 2 /|,r 2 + . . ., 
the second derivative of the ACF needed to compute the LCR 
in ([H) is evaluated as: 



p Ray (0) 



-4-7T 



EN j2 fl 
i=l 1 i J D 



(13) 



Hence the three parameters, r, 9 and i?(0), are available and 
^ gives the approximate LCR. 

B. LCRs for Rician Fading 

As in the Rayleigh fading case, the instantaneous aggregate 
interference, I R i C {t), for this scenario is given as: 



TV 



i R ,c(t) = Y, I i\Ht)f 



(14) 



where hi (t) is Rician and N, I\ , I2, ■ ■ ■ , In are as defined 
in ©. Hence, I R i C (t) is a weighted sum of noncentral chi- 
square (% 2 ) random variables. Using the same approximation 
philosophy as that used in the Rayleigh case, we propose 
approximating ( TT4T > by a single non-central \ 2 ■ This approach 
is less well documented but has appeared in the literature (see 
[13]). Also note that a scaled, rather than a standard, non- 
central x 2 distribution is required for fitting. A noncentral % 2 
variable with v degrees of freedom, non-centrality parameter 
A and scale parameter a has the following PDF: 



P( x ) = o ex P 



-(* + «*) V«?f) - 7_(VA^), (15) 



where I/ v -2)/2 i s a modified Bessel function of the first 
kind with order (v - 2)/2. Fitting the PDF in (Q3]i to the 
variable in ( TPfl i is performed using the method of moments 
technique so that the approximate noncentral % 2 has the 
same first three moments as I R i C (t). Note that there can be 
numerical difficulties with the approach for certain values of 
7i, I2, ■ ■ ■ , In- Results are shown in Sec. IV for cases where 
the estimation procedure was successful. Further research is 
necessary to make the methodology robust to all possible 
interference values. 

Next we consider the LCR of the noncentral \ 2 process 
which is used to model I R i C (t). Consider a generic scaled 
noncentral \ 2 process, g, given by: 

g= SkCj.il*)., (16) 

a 

where Y^i=i = A is the non-centrality parameter, v is the 
order (degrees of freedom) and a is the scale parameter. Note 



that for simplicity we have omitted the dependence on time so 
that g(t) is denoted by g. The variables Xi, X2, ■ ■ ■ , X v are 
i.i.d. Af(0, 1). Using the basic formula of Rice, the LCR of g 
across a threshold, T, is given by: 



1 

gPg\g(g\T)p g (T)dg, 



(17) 



where p g (.), Pg,g(-,-) and P g \ g (-) are the PDF of g, the 
joint PDF of g and g and the conditional PDF of g given g 
respectively. The derivative of ( [ToT l with respect to time gives: 



El=i 2(X 4 + 6i)Xi 



a 



(18) 



Now, using ( TTol l and HSi . p g \ g (g\T) is a Gaussian PDF corre- 
sponding to the distribution jV*(0, 4-fT/a) with 7 = Var(X). 
Thus, the LCR in (T7\ becomes: 



Pg( T )Va 

V^nT Jo 

Pg(T) 



gcxp 



-gg 

87T 



dg 



2 7 T 



(19) 



Assuming the Jakes' fading process for the Gaussian 
variables, X\, X2, ■ ■ ■ , Xm, each Xi has ACF given by 
J (2irf D T) = 1 - tt 2 / 2 r 2 + .... Hence, 7 = Var(X) -> 
27T 2 /!, as t ^ [14]. Substituting, p ff (T) from (B), the 
LCR becomes, 



LCR; 



Pg(T) 



.(VAoT). 

(20) 



Note that this is the LCR of the process in (fTol l. and appears 
to be a new result. However, the noncentral x 2 which fits 
Imc(t) will almost certainly not have an integer order. Hence, 
this derivation for integer ordered noncentral x 2 is applied to 
the fractional order case of interest. At present the validity of 
this approach is only a conjecture, but the results in Sec. IV are 
very encouraging. Note that a similar extension for a central 
X 2 with integer order [15] to a central x 2 w ith fractional 
order [12] has been shown to be correct. A comparison of 
the theoretical results given in ^ and ( f20b with Monte Carlo 
simulations along with a commentary on the results is given 
in the next section. 

IV. Results 

Throughout the section we assume the following parameter 
values: shadow fading variance, a = 8 dB, path loss exponent, 
7 = 3.5, radius of PU coverage area, R = 1000 m, radius of 
CR coverage area, R c = 100 m, CR density, 1000 CRs per 
square kilometer, an activity factor of 0.1 and fr> = 25 Hz. 

A. Imperfections in the REM 

In practice, the radio environment is often modeled by 
dividing an area into a regular grid (typically composed of 
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Fig. 1 . Interference CDF for an REM enabled CR network for several values 
of A and decorrelation distance, Dd = 100 m. 
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Fig. 2. REM grid size, A, vs decorrelation distance, D^, for different 
threshold buffer sizes. 



100 m x 100 m grid boxes) and assuming that the fading 
conditions in any grid box can be approximated by a single 
point at the center of the box. For example, drive testing of 
cellular networks to validate path loss models and predicted 
signal coverage follows this approach. Clearly, larger grid 
sizes result in errors between measurement and prediction. On 
the other hand, reducing the grid size results in a large data 
overhead. Figure 1 shows the cumulative distribution function 
(CDF) of the magnitude of the actual CR-PU interference 
when the REM is estimated via a grid size ranging from 1 
m x 1 m to 100 m x 100 m. The REM approach aims to 
maintain a 2 dB SINR buffer for the primary, but this is only 
possible with a perfect REM. When A = 1 m the 2 dB buffer 
is nearly achieved but for a grid size of 50 m x 50 m, the 
interference exceeds 3 dB for approximately 8% of the time. 
For a grid size of 100 m x 100 m, 3 dB is exceeded 30% 
of the time. In effect this means that if REM information is 
derived from a coarse grid, the buffer size must be increased 
or the CRs must back off from the buffer. 

The effects of increasing the buffer or backing off the CRs 
are shown in Figs. 2 and 3 respectively. In Fig. 2 the PU has 
a target 2 dB buffer but due to the imperfect REM it will not 
always be achieved. Hence an extra buffer is permitted beyond 
which the CRs are only allowed, 5% of the time. In Fig. 2 





Fig. 4. LCR results for different fading conditions. The solid lines represent 
analytical results. Simulation values are shown by the circle, star and triangle 
symbols. 



Fig. 6. LCR results for the dominant and no dominant interferer cases in 
a Rician (K = 10 dB) fading scenario. The solid lines represent analytical 
results. Simulation values are shown by the circle and star symbols. 



this scenario is denoted by the legend, Threshold = (original 
+ extra) dB. The effects of spatially correlated shadow fading 
are also considered in Fig. 2. Shadow fading is correlated over 
any given area and the level of this correlation has a simple 
effect on the REM grid size. For highly correlated areas a 
coarse grid (large A) will be acceptable whereas in areas of 
low correlation, a fine grid (small A) will be required. Figure 
2 shows the REM grid size vs the decorrelation distance of 
the shadow fading. For a given interference degradation (say 
the buffer value plus an additional 2 dB) a large decorrelation 
distance (say 500 m) enables a coarser grid size 165 m x 165 
m relative to a decorrelation distance of 100 m (typical for 
dense urban areas) when the grid size is 70 m x 70 m. 

In Fig. 3 we consider a back off in the CR allocation policy. 
In order to meet the nominal 2 dB SINR buffer at least 99% 
of the time, the CRs have to target a reduced buffer which is 
less than 2 dB. Figure 3 shows this buffer vs A for various 
values of D^. For a grid size of 25 x 25 m and a decorrelation 
distance of 100 m, the interference buffer is 1.05 dB. Figures 2 
and 3 are instructive in determining the grid sizes for different 
radio environments. 



B. LCR and AED of CR-PU Interference 

Figures 4, 5 and 6 show the LCR (normalized by Doppler 
frequency) of the interference for different types of fading and 
interference profiles. The interference profiles, i.e., the values 
of Ji, I2, . . . , In> are determined for each fading type via a 
simulation of the CR allocation policy. From 1000 simulations, 
two sets of interferers are selected. The first set has a dominant 
interferer and corresponds to the set of interferers with the 
highest variance. The second set has no dominant interferer 
and corresponds to the set with the least variance. 

For all types of fading, the maximum LCR is observed close 
to the buffer value. This is because the CR allocation method 
gives a mean interference level close to the buffer. Even in 
strong LOS conditions (K = 10 dB), the interference shows a 
significant number of level crossings across the buffer due to 
the scattered component. Figure 5 shows the case of Rayleigh 
fading where the interference budget is dominated by a single 
large interferer with a number of smaller additional interferers. 
Also shown is the case where no dominant interferer exists. 
Figure 6 shows the same results for a Ricean channel with 
K = 10 dB. Figures 5 and 6 show that when there are 
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Fig. 7. AED results for the dominant and no dominant interferer cases in 
a Rician (K = 10 dB) fading scenario. The solid lines represent analytical 
results. Simulation values are shown by the circle and star symbols. 

many small interferers, the resulting interference is more stable 
compared to the dominant interferer case. Furthermore, the 
Ricean channel always returns a lower LCR as compared to 
a Rayleigh channel. The results in Fig. 6 are quite promising. 
Under the near LOS conditions that may be present with 
small cell radii, the CR-PU interference has a much lower 
level crossing rate across the interference buffer for the no 
dominant interferer case. Hence it may be a desirable part of 
the CR allocation policy to avoid any single user which takes 
up a significant part of the buffer. Finally, for completeness, 
we show the AED results corresponding to Fig. 6. The 
AED follows from the LCR using standard results [16]. As 
expected the time spent by the interference above a threshold 
decreases as the threshold value increases. Therefore, for the 
no dominant interferer case, the interference seldom crosses 
the threshold (see Fig. 5), and when it does, it only exceeds 
the threshold for a small period of time. Finally, we note that 
all figures show an excellent agreement between the analytical 
approximations and the simulations. 

From the point of view of PU system designers, the follow- 
ing questions are important: 

• How much is the CR-PU interference? 

• Can it be controlled? 

• How often will it exceed a threshold? 

• What happens when it does exceed the threshold? Does 
it stay above for a long time or quickly return back to 
acceptable levels? 

• How do all of the above change with the type of fading? 
Figures 4-7 shed an interesting perspective on all the above 
questions. 

V. Conclusion 

In this paper we have shown that interference degradation to 
the PU can be significantly underestimated if the channel state 
information needed to estimate interference levels is derived 
from a coarse REM. For practical deployments, this may mean 
that the PU has to accept a much larger interference from 
the CRs or the CRs may need to set a more conservative 
interference target. This will reduce the number of CRs 



allowed. We also determine the LCR and AED for the CR-PU 
interference and show that the maximum LCR occurs close to 
the maximum allowed interference level for both Rayleigh and 
Rician channels. The LCR results show that it is desirable for 
the interference to be made up of several small interfering CRs 
rather than a dominant source of interference. The LCR of the 
former case is more stable than the latter. The AED results 
also show that the interference exceed the threshold value for 
small periods of time in the latter case. 
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